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Abstract 

We consider the M/M/l queue with processor sharing. We study the 
conditional sojourn time distribution, conditioned on the customer's 
service requirement, in various asymptotic limits. These include large 
time and/or large service request, and heavy traffic, where the arrival 
rate is only slightly less than the service rate. The asymptotic formulas 
relate to, and extend, some results of Morrison [4] and Flatto [8]. 
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1 Introduction 



One of the most interesting service disciplines within queueing theory is pro- 
cessor sharing (PS). Unlike other disciplines such as first-in first-out (FIFO), 
under the PS discipline every customer gets a share of the server. The ad- 
vantage of this discipline over FIFO is that customers requiring only short 
amounts of service get through the system more rapidly than with other ser- 
vice disciplines. The PS discipline was apparently introduced by Kleinrock 
PQ, [2], and has been the subject of much further investigation over the past 
forty years. 

In [3] , Coffman, Muntz, and Trotter derived an expression for the Laplace 
transform of the distribution of the waiting time in the M/M/l-PS model. 
We shall denote the waiting time (in the steady state) by W and the sojourn 
time by V, which is the total time from when a customer arrives to when that 
customer leaves the system, after being served. In the M/M/l-PS model we 
denote the Poisson arrival rate by A and the exponential service time density 
by pe~^ x . The traffic intensity is p = A//i. If we are given a service time x, 
then we let V(x) be the conditional sojourn time. The "waiting time" W(x) 
is defined by W(x) = V(x) — x, and with this notation E[e~ sW ^- ) ] is the 
Laplace transform of the conditional waiting time distribution derived in [3] . 

Using the results in [3], Morrison [1] studied the unconditional sojourn 
time distribution in the M/M/l-PS model, in the heavy traffic limit where 
A t P (thus p | 1). Let us denote the unconditional sojourn time density by 
p{t)dt = Pr [V G (t, t + (it)] and the conditional density by 



and we note that p(t\x) has support for t > x, with a probability mass along 
t = x. 

In jl] Morrison assumes that A f p and defines e = 1 — X/ p = 1 — p. 
He obtains approximations to p(t) on the time scales t = 0(1), t = 0(e _1 ), 
t = 0(e~ 2 ) and t = 0(e~ 3 ). Most of the mass, in this asymptotic limit, is 
concentrated in the range t = 0(e _1 ). 



p(t\x)dt = Pr [V(x) e(t,t+ dt)] . 



(1.1) 



For the M/M/l-PS model we have 




(1.2) 
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A service discipline seemingly unrelated to PS is random order service 
(ROS), where customers are chosen for service at random. The M/M/l- 
ROS model has been studied by many authors, see Vaulot [5], Riordan [6], 
Kingman [7J and Flatto [8] . In [8] an explicit integral representation is derived 
for the waiting time distribution, from which the following tail behavior is 
computed: 

Pr [W ROS > t] ~ e- At e- Btl/3 Cr 5/e , t -> oo. (1.3) 

Here A = (1 — ^fp) 2 ^ we scale time to make the service rate p = 1, and 
p < 1. This formula appeared previously in the book of Riordan [9] (pg.105) 
and was apparently first obtained by Pollaczek [10] . Cohen [11] established a 
relationship between the sojourn time in the PS model and the waiting time 
in the ROS model, 

Pr[V>t] = ^Pr[W ROS >t], (1.4) 

which extends also to the more general G/M/l case. In [12] relations of 
the form (11.41) are explored for other models, such as finite capacity queues, 
repairman problems and networks. Later in this paper we discuss the rela- 
tionship of the tail formula in (I1.3P to the heavy traffic expansions in [I]. 

In this paper we study the conditional sojourn time distribution for the 
M/M/l-PS model in two cases. First we consider a fixed p < 1 and obtain 
expansions of p{t\x) for t and/or x — > oo. From these (11.30 is readily obtained 
by using (II. 2p for t large. Then we consider the heavy traffic limit where p f 1, 
and again obtain approximations for several ranges of the space-time plane. 
From these results all of the expansion in [I] can be recovered. 

We mention some other work on PS queues. The G/M/l-PS model was 
studied by Ramaswami [13] and the M/G/l-PS model by Yashkov [14J, 
[T5] and by Ott |16j . For the latter there is a complicated expression for 
the Laplace transform of the conditional sojourn time distribution. For 
the special case of M/D /1-PS this expression simplifies considerably, and 
then the tail behavior was obtained by Egorova, et. al. [17] in the form 
Pr[V > t] ~ B'e~ 1 . Comparing this to (11.41) with (11.31) the interesting ques- 
tion arises as to what are the variety of possible tail behaviors for the general 
M/G/l-PS model. For the G/G/l-PS queue the tail was shown by Mandjes 
and Zwart [18] to be roughly exponential, in that log{Pr[V > t]} ~ — A$t as 
t — ► oo. This assumes that the arrival and service distributions have expo- 
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nentially small tails, so that their moment generating functions are analytic 
in some part of the left half-plane. 

There have also been several recent investigations into tail behaviors of 
PS models that are "heavy tailed", in that the service time distribution has 
an algebraic tail. The M/G/l model with this assumption is analyzed by 
Zwart and Boxma [19], where it is shown that the sojourn time has a similar 
algebraic tail. A thorough recent survey of sojourn time asymptotics in PS 
models can be found in [20] , where both heavy and light tailed distributions 
are discussed, as well as various approaches to obtaining the asymptotics. 

The remainder of the paper is organized as follows. In section 2 we 
summarize the main results, as Theorems 1 and 2. In section 3 and 4 we 
sketch the main points of the derivations. A brief discussion appears in 
section 5. 



2 Summary of results 

We denote the mean rate for the Poisson arrival process by A, and assume 
that the mean required service time is unity. Then the traffic intensity is 
p = A/p = A > 0, where p = 1 as in the notation of Coffman, et.al. [3]. 

It was shown in [3], under the stability condition p < 1, that the Laplace 
transform of the equilibrium waiting time distribution, conditioned on the 
job requiring x units of service, is given by 



-p(l— r)x 



E r p --w(x)i = (1 -p)(l -pr 2 )e- 

1 J (1 -pr) 2 -p(l _ r )2 e -x(l-pr2)/r ' 



(2.1] 



where 

r = r(s; p) = — \(1 + p + s) - + p + s) 2 - 4pj . (2.2) 

Ip L J 

Taking the inverse Laplace transform of ( 12. 1ft . the probability density of the 
sojourn time, conditioned on service time x, is 

p(t\x) = — [ e s{t - x) E[e- sW{x) }ds 
2ni J Br 

1 C e st e -(*+p-pr)*(i -p)(l - pr 2 ) 
2m J Br (1 - pr) 2 - pjl - r )2 e -(i-^ 2 )A (2 ' 3) 
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Here Br is a vertical contour in the complex s-plane, on which 3?(s) > 0. Note 
that this density contains the term 8(t — x)(l — p)e~ px , which corresponds 
to probability mass at t — x. The sojourn time is equal to the service time 
if the customer arrives to an empty system and there are no arrivals during 
that customer's service time. Analyzing the integral f)2.3p . we obtain the 
following results: 

Theorem 2.1 For p < 1, the conditional sojourn time density has the fol- 
lowing asymptotic expansions: 

1. x — > oo, t — x — > + with x(t — x) = 0(1), 

p(t\x) ~ (1 - p)JJ^-e-' at I 1 (2y/px(t-x)) (2.4) 

where I±(-) is the modified Bessel function. 

2. x,t — > oo, 1 < t/x < oo, 



(1 - pW/4 e 2Vp(l-*/*) e -(l+/>V , , 

p(t\x) ^ ^ — (2.5) 

2^^/x [(t/x) (t/x - 1)] 3/4 y/t/x ~ y/p(t/x - 1) 

3. x -> oo ; t/x 2 = a = 0(1), 

(1 + ^p) e -^-^ H e^-^ x 



p(t\x) 



20Fp 1/4 (l - v/p)a 5 / 2 x 3 

oo 

• e-( 2n+1 ) 2 vW(4a) [ (2n + if^p _ 2a] . (2.6) 

n=0 

By the Poisson summation formula, we also have 

(2.7) 

4- x — 0(1), t — > oo ; 

p(t|z) ~ F(x)e r * (x)t , (2.8) 
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where for x > = 

. , 2./p(l - p)e^-^ cosv ^sm 2 vcos( y ^/pxsmv) 

Fix) — 

(1 + p — 2y / pcosf ) 2 v /pa; cosf + (1 — p)(l + p — 2y / pcost>) 

• [(1 — p) sinvtan(y / pa;sinf ) + 2 v /p — (1 + p) cosi>] 

and 

r*(x) = -1 -p + 2^/pcos(v(x)) (2.9) 
where v — v(x) is the smallest positive root of 

2 



(2.10) 



For < x < we nave 

r*(x) = — 1 — p — 2 v /pcosh(-u(a;)) 
where u(x) is the unique positive root of 



n+^pe- u \ 
V Vp + e ~ u ) 



2 

_ 2y/~pxsinhu 



and in F(x) we must replace (sin v, cos v) by (isinhu, — cosh -a). When 
x = x* we have r*(x*) = —(1 + ^fp) 2 , v(x*) = n, u(x*) = and 

6ygl + y^) 2 / 1 

— 1 : — zr^ CXP — — — 1 



r*(x) aas £ne following expansions, as x — > oo: 



and as rr — > + : 

r*(x) = 



logp 



(2.12) 



2p(logp) 2 + (p 2 -l) logp- 4(1 -p) 2 

+ (logp) 3 X + C^j. 
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In the asymptotic matching region between cases 3 and 4 we have 
2tt 2 (1 + y/p) 



p(t\x 



•exp <^ (1 - y/p)x 



(2.13) 



: i -v / p)-^ + 



vr 2 , 2tt 2 (1 + Vp) 



^Tir 2 p(l - y^x 3 

which holds for t — ► oo and x — > oo, with x = 0(t 1//3 ) or t 1 / 3 C i < i 1//2 . 
By removing the condition on x and noticing the relationship (jl.4p between 
processor sharing and service in random order, we can recover Theorem 1.2 
in Flatto [8]. 

For the heavy-traffic case, where p is close to 1, we let e = 1 — p — > + . 
Then we have the following results: 

Theorem 2.2 For p = 1 — e, where e — > + , we let t = T/e and x = X/e. 
The conditional sojourn time density has the following asymptotic expan- 
sions: 



1. x = 0(1), t = 0(1), 

p(t\x) ~ 

where 

2. x = 0(1), T = 0(1), 

p(t\x 



e st e ~(s+l-r )x ^ _ r j2^ 

2m J Br (1 - r ) 2 (l - e- x ^- r o^ r °) 



ds, 



(2.14) 



2 



s + 2 - Vs 2 + As 



— e - T / x — x ^ 
a; 6 



2x — T 



-T/x 



e 2 + 0(e 3 ). (2.15) 



5. X = 0(1), T = 0(1), T - X -> 0+ mft (T — X) = e 2 ^ = 0(e 2 ), 

p(t|rr) ~ e" x/e e x sfxjTJx (2 ) . (2.16) 
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I X = 0(1), T = 0(1) and 1 < T/X < oo, 

e 3/2 (VT+v/T^X^ 2 

p(t\x) 



cxp 



2v^FT 3 / 4 (T - X) 3 / 4 
expji 2^T(T - X)+ X - 2T j 



T - y/T(T - X) 



(2.17) 



5. X = = O(v^), T = 0(1), 

2e 3/2 

p(*f) r 



'txT 



n=0 



3 -(2n+l) 2 Z 2 /(4T) 



(2.18) 



and fry £/ie Poisson summation formula, 

e 3/2 

n=l 



p(t|x) 



[l + 2 ^(-l)« e -™ 2 ^]. 



(2.19) 



5. X = O(l) ; T = 0/e = O(e 1 ) ; we owe £/ie expansion in three different 
forms: 



(a) 



p(t\x) 



4e 2 e x/2 e -e/4 



ni J Br (1 + 2^) 2 - (1 - 2v^) 2 e- 2 ^ 



2n 



p(t|z) 



(2n)! 



n=0 



Z=0 



/!(2n-/)! 



D-i(z n ) - AD^i_i(z n ) + V2BD_ l _ 2 (z n ) 



(2.20) 



(29) i/2 



(2.21) 



ifere D v (-) is the parabolic cylinder function, and z n = ( - 2n+ ^ / ^ +e - 



(c) 



p(t\x) ~e 2 J2 F i X ; v n )e r * (x > v " )e . (2.22) 



n=l 



Here 




\2iv - 1 ) 




(2.23) 



If Q = cr/e = 0(e x ) ; we have 



p(t\x) ~e 2 F(X;i; 1 )e r - (X:uiW£ 



(2.24) 



where V\ = Vi(X) is the unique root in the interval (0, ir/X) of A2.23\) . 

By removing the condition on x and using the results ( 12. 24ft and (12.15)1 . 
we can recover the results in Morrison [4]. The details are given in section 
4. Here we also relate the unconditional tail expansion of Flatto [8] to the 
heavy traffic result of Morrison [4] . 

3 Brief derivations for the case p < 1 

We assume that the traffic intensity p is fixed and less than one. Con- 
sider the function (12.11) with (I2.2p . We first observe that if we replace 
- a/(1 + P + s) 2 - 4p by +y/(l + p + s) 2 -4p in (Q, then r becomes l/(pr). 
But ( 12. ip is invariant under the transformation r — > l/(pr). From (12.21) and 
the fact that r(s; p) appears in (12.11) it would seem that the integrand in (12.31) 
has branch points at s = s± with = — (1± ^fp) 2 = — 1 — p=F2 v /p, but from 
our discussion it follows that (12. ip has in fact no branch points. It has poles 
where the denominator in (12. ip vanishes. We also note that r(s_; p) = 1/v/p 
and (1 — pr 2 )/r > for s real and s > s_. 

We define 




(3.1) 
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and rewrite (12. 3p as 

(l-p)e-^ / a „ w 1-pr 2 

WJ 27TZ A/ (l- pr )2_ p (l_ r )2 e -x(l-p r2 )/r CtS 

~ (1 = ^ I e^ l - pr2 ^ P \ds (3.2) 
2« 7 Br [l-pr(s;p)] 2 

where the last approximation holds when x — > 00 and 5J(s) > s_ on the 
Bromwich contour I?r, which is a vertical contour in the complex s-plane. 

Consider first the limit x,t — > 00 with 1 < t/x < 00. The second inte- 
grand in (13. 2p has saddle points where <p'( s ) = so that there is a saddle 
along the real axis at s , where 

so = s - = -1 - p+ (3.3) 

Note that s_ < s < 00 whenever 1 < t/x < 00, and s — > 00 as t/x J, 1, 
while So — * s_ as t/x — > 00. If we shift Br in (13. 2p to Br', on which 3?(s) = sq, 
and use the standard Laplace method (see, e.g., Wong [21]) we get 

pm „ (l-,)e-~[l-^(,. ; , )]«■*«> 
V27T [1 - pr(s ;p)] a/x<p"(s ) 

But from (jHHJ) and (Q we get 



r ( s o;p) = x/ 1 —^, (3-5) 



x<p(s ) = (a: - t)(l + p) + 2y/t(t - x), (3.6) 

and 



*>o) = ^,3 • (3.7) 



2t 3 / 2 (t - x ) 3 / 2 
^x 3 

Using ([SSD-dXTD in ([S3D leads to fl2~3|) . 

Next we consider x, t — > 00 but with t/x ~ 1. The previous calculation 
is not valid since the saddle point so — > 00. Now we return to the second 
expression in (13.21) and note that r(s) ~ 1/s as s — > 00. Then we approximate 
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the integrand for s large (more precisely we can scale s = 0(x) with x(t—x) 
0(1)) to get 



p(t\x) 



:i-p)e- 

2m 



px 



At-x) e px/s ds 



(3.1 



Br 



1 - p)e- px 5{t -x) + 



[1 - p)e~ px 



Br 



[l-p) e -P*S(t-x) + (l-p)e-'"-£ 



s{t-x) ^px/s _ ds 
e s(t-x) 



1 

2vri 



Br 



( 



jpx/s 



1) ds 



The last integral is equal to J ^2a/px(£ — x)J — 1, where J is the modified 

Bessel function. Differentiating this with respect to t leads to (12.41) for t > x. 

Now we consider x, t — > oo but with x/t small. In this limit the sad- 
dle point So — > s_, which is a branch point of the integrand in the sec- 
ond expression in (13. 2p (but not of the first). We shift the Br contour 
to 5r_, which lies slightly to the right of s = s_. As s j s_ we have 
x(l — pr 2 )/r ~ 2p 1//4 v /s — s_x, so for x large we scale s — s* = 0(x~ 2 ) by 
setting 

s = -I- p + 2^p + z/x 2 = + z/x 2 . (3.9) 

We then have 



xip(s) = s(t — x) + p 



x 



P 



-3/4 



^~z + 0(x- 1 ) 



(3.10) 



= -(1 - y/pft + (1 + p - y/p)x - p 1/4 v^ +~l z + Oix- 1 ). 

If we assume that x = 0(\rt) with a = t/x 2 , then ( 13.21) becomes 

p(t|z)~ 2p~ 1 / 4 (l - p) e -(Wd a t e (W*)* (3.11 

-dz. 



27TZ 7 Br (1 - ^)2 [i _ e -2p^] 
Expanding ^1 — e~ 2pl/4v/ ^ as a geometric series and using the identity 

(3.12) 



1 

2tH 



Br- 



cxp 



(2n+l) 2 > 



4a 



[(2n + l) 2 v/p - 2a] a" 5/2 
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we obtain (12. 6p . Then (12.71) follows upon using the Poisson summation for- 
mula in the form 



oo oo oo /»oo 

E = E ^ 27rm ) = E / e^VWl/, (3.13) 

n=— oo m=— oo m=— oo ~ 00 

where ^ is the Fourier transform of and 7/>(n) is given by the right side 
of (I3.12j) . Alternately, we can set z = w 2 in (13.111) and after some contour 
deformation write the integral as 



w 2 g a«; 2 e -p 1 / 4 w) 



7ri(l - ^p) 2 x 3 J Brjr 1 - e - 2 P 1/4 



dw, (3.14) 



where dt(w) > on Br*. The integrand in (13. 14j) has simple poles at w = 
p~ 1 / 4 N7ii for N = ±1, ±2, and is anti-symmetric under the transformation 
w —>■ —w. We can thus write it as an infinite residue series and this again 
leads to (TO) . 

Finally, we consider the CclSG X = 0(1) and t — > oo. We use the conformal 
map s = —1 — p + 2 v /pcosh?y. Then in the 77-plane, f !3.2j) becomes 

Mx)= h^zA (3.15) 

Z7TI 

e (-l-p+2 v / Pcoshr?)t e (l-2 v ^coshr?+ v ^e- , ')xQ _ e ~2v\ gi n h 77 

— drj. 



c 



Here the contour C goes from 77 = 00 — 77r/2 to 77 = 00 + z7r/2 as indicated 
in Figure 1. It cannot enclose any of the poles, whose location is discussed 
below. We let 77 = u + iv where u = ^(77) and v = ^s(rj). 

Note that 77 = is a simple zero of the integrand, since (1 — e~ 2r] ) sinh 77 = 
0(t] 2 ) and the denominator in ( 13.151) has a simple zero at 77 = 0. From nu- 
merical calculations we find that the dominant singularity (which maximizes 
— 1— p+y/pcosliT?) occurs either along the imaginary axis for |i>| < tt, or along 
the lines 3 (77) = v = ±7r. In Figure 2 we sketch all the singularities of the 
integrand in (13 . 151) for x — 1, p — 0.3 (so that x > x* = ^i^^ = 0.5335). 
The dominant singularity is now on the imaginary axis at 77 = ±iv where 
v = 1.6111. In Figure 2, the singularities are the crossing points of two 
curves, corresponding to the real and imaginary parts of the denominator in 
(13.151) . with 77 = u + iv. As x decreases through the value x*, the dominant 
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singularity "turns the corner" and begins to lie on v = ±7r. In Figure 3 we 
sketch all the singularities of the integrand in (I3.15P for x = 0.5, p = 0.3 (so 
that x < x* = 0.5335), with the numerical value of the dominant singularity 
being rj = ±in+u where u = 0.5195. In Table 1 we give the dominant singu- 
larity (with v > 0) for various values of x, and various p < 1. When x = x*, 
the dominant singularity is at rj = ±in. Note that when x = x*, rj = in is a 
double zero of the numerator and a triple zero of the denominator of (I3.15P , 
so f] = in is a simple pole of the integrand. The Table also indicates the 
behavior of the singularity for x — > and x — > oo, where (12.101) and (12.111) 
apply. 

Expression (12. 8p follows by indenting the contour in Figure 1 around 
the two dominant singularities and using the symmetry (v — > —v) of the 
integrand in (I3.15p . To obtain (12.111) we use (I2.10p and note that v — > as 
x — > oo with Xy/psinv — > n, so that v ~ n/(y/px), which can be refined to 
the estimate 

tt 7r ( 1 -I- . f~n\ 

0(x' 3 ). (3.16) 



Vp x Pi 1 ~ \fp) %2 

To obtain (12.111) we use 

s = -1 - p + 2 v ^cosw = -(1 - ^) 2 - ^/pV + 0(v 4 ) 

and (13TT6D . 

To get the dominant singularity as x —>■ + we set rj = u + 27r to get 



Vp + e- u 



(3.17) 



As x — > + we must have w — > oo and, using 2sinhw ~ e u , we obtain from 

(TO 

Vpa;e u = log[l/p+0(e- u )] (3.18) 
so that ^ 

M ~ - log X + log log(p _1 ) - - log p. 

By considering higher order terms in (13.181) this can be refined to 

u= -logx + loglog(p _1 ) - ^logp - — ^x 

2 (logp) 2 

, p(logp) 2 + (p 2 - 1) logp - 6(1 -p) 2 ^ 2 + 3 



(logp) z 
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Then ( 12.121) follows by using r*(x) = — 1 — p — ^fp{e u + e "). 

To compute the unconditional density p(t) we use (11.21) with \i = 1. The 
major contribution will come from the asymptotic matching region between 
the scales x = 0(1) and x = 0(\fi). In this range we use (12. 8p with r*(x) 
approximated by (12.111) and 



F(x) 



p(l - ^p)x 3 



(i-yp)z 



x 



oo, 



(3.19) 



which follows from the definition of F below (12.81) and from (I3.16p . Scaling 
x = t l l 3 y = Oit 1 / 3 ), (ll.2p becomes asymptotically 



7[- 



-2/3 



t 2/3 

^ot- 1 / 3 y 



i-Vp)p 

" 2vr 2 (l + yp) 

A 1 - Vp)v 3 . 



[3.20) 



1 

-3 eX P 



exp 



t 1/3 -V~pv 



7T 



dy. 



Here Xq > so as to avoid integrating through y = 0, and Xq 3> 1. 

Expanding (I3.20p by the Laplace method, with the major contribution 
coming from y = {2tx 2 / p) 1 ^ 3 , leads to 



p(t) ~ e- At e' Btl/:i C^ 6 



(3.21) 



where 



.4 



;i - VP) 2 , 



a = 2 2 / 3 3- 1 /vv 5/ 

By using f[3~2Tj) and ([13]) we get 

p(*) = -^Pr[V>*p 



12 



VP 



exp 



~ l AC*r 5 ' & e~ At e- Btiri 



and this agrees with the result of Flatto (see Theorem 1.2 in [8]), if we allow 
for the different scaling of time used therein (where the arrival rate was taken 
as unity). 
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4 Brief derivations for the case p « 1 

Now we consider the case in which the traffic intensity is close to one, and 
let p = 1 - e with < e < 1. 

First, we consider x = 0(1) and t = 0(1). From (12. 2p we have 



This leads to (12.141) . On this scale the solution does not simplify much, but 
there is little probability mass in heavy traffic on the time scale t = 0(1). 

Next, we consider x = 0(1) but large time scales t = T/e = 0(e _1 ). In 
(12.31) . we replace p as 1 — e and scale s as ew, and we have 



Note that by removing the condition on x and expressing the result in terms 
of the modified Bessel function, we can recover the result of Morrison (see 
(3.12) in [1]). The term proportional to S(T) does not actually mean there is 
mass at T = 0, but rather corresponds to the small (0(e)) mass that exists 
in the shorter time scale t. 

Now consider x = X/e = 0(e^) and t = T/e = 0(e^) with 1 < T/X < 
oo. We define 





Then (12.151) follows upon using the following identities 





4(8) 



T 



l) + (l-e) (r(s-l 



X 




Ms) + Ms)e + 0(e 2 ) 
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and rewrite (12. 3p as 

e f e^ x ' e [1 - (1 - e)r 2 ] 
n(t\x) = / - - - - ds 

V ' ; 2m J Br [1 - (1 - e)r] 2 - (1 - e)(l - r )2 e -a.[l-(l- e )ra]/r 



(1 -e)r 



2 



re 



*(")*A ds 



2m J Br [1 - (1 - e)r] 

„ / .f^A e ^ x e^ x ^da. (4.1) 

27riJ Br \ s 1 ; 

The saddle point in the last approximation in (14. ip . where 4>' (s) = 0, is 
s = Sq where 

y 2T X 

Note that in the heavy traffic case, s_ = —(1 — ^fp) 2 = 0(e 2 ), s_ < Sq < oo 
whenever 1 < T/X < oo, Sq — > oo as T/X | 1, while Sq — > as T/X — ► oo. 
By shifting _Br in ( 14.11) to Br", on which 5J(s) = Sq, and using the standard 
Laplace method, we have 

p M x ) ~ ^V ^ + V^ >K)x >(s g)x/ e 

This leads to (12TTTD . We note that if X -> oo and T -> oo but T/X = 0(1), 
the approximation (I2.17P remains valid. 

For the case X = 0(1), T = 0(1) and T — X — > + , we notice that the 
saddle point Sq — > oo. Then by the same argument as in the case when p is 
fixed and less than one, we can easily get (12.161) . where T* = (T — X)/e 2 = 
0(1). 

Next, we consider X = = 0(- v /e), T = O(l). It follows that T/Z 2 = 
0(1) and the saddle point Sq — ► 0. We shift the contour Br to -Br , which 
lies slightly to the right of s = 0. By scaling s = 0(e) = ew in the first 
equation in (14. ip . we have 

2e 3 / 2 [ e wT e~^ z 

p(t\x) ~ / — -. „ _ . aw. 

FK 1 ' 2ni J Bro yfij](l- e ~ 2 ^ z ) 

Expanding (l — e~ 2 ^" z ^ 1 as a geometric series and using the identity 

1 r c wT c -(2n+\)JwZ 1 

1 ' e e = 1 c -(2n+l)^V(4T) 
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we obtain (12.181) . Then (12.191) follows by using the Poisson summation for- 
mula (13.131) . Note that from (12.191) . we can easily verify that cases 5 and 2 
in Theorem 2 asymptotically match, in the intermediate limit where x — > oo 
and Z — > 0. Similarly, cases 4 and 5 match in the limit X — > and Z — > oo, 
which follows easily from (12.181) . 

Now we consider X = 0(1) and T = O/e = 0(e _1 ) (thus x = 0(e~ x ) and 
£ = 0(e~ 2 )). In the first expression in (14. ip . we scale s = 0(e 2 ) by setting 
s = (£ — l/4)e 2 , which leads to (I2.20p after some calculation. 

Furthermore, we expand the integrand in ( 12.201) as a geometric series, and 
we have 

1 1 ; ™ v.(i + 2v^) 2 ^ vi + 2ve; 

(4-3) 

Note that if we let X — > oo, then the n = term in H4.3[) dominates, and we 
have 

2«y B? , (i + 2^) 2 ' 



e- /«*> y| e _ x2/(40) _ 4 + X + 9 e6/Wfc /X + 9 



Here we used 



4\/7re Wn 8 V 2V6 

-X 2 /(40) 

X — ► oo. 



2 /" 

erfc(^) = 1 — erf(z) = —= / 



oo 

-t 2 Jj. L -z 2 

e at ~ e , as 2 — > oo. 

V7TZ 



Then ( 14 . 3 1) becomes, for fixed and X — > oo, 

p( f | x ) „ JfL e x/2 e -e/4 e -xV(4e)_ (4 _ 4) 
V7r6 

When X = 0(e- 1 ) and T = 0(e _1 ) but with T/X = 0(1), (l2~T7l) remains 
valid, and letting T/X — > oo in ( 12.171) regains ( 14.41) . This again verifies that 
these two cases asymptotically match. 

We return to ( 14. 3p . and let \/£ = (z — l)/2, with which the integral 
becomes 

p 0/4 (■ 00 / n\2n 

^_ / ( z _ 1)2 y \± J± e (2n + X)X/2 e -K2n+l)X+eW2 e Z V4 d ^ 
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Here the contour C can be taken as the imaginary axis in the z-plane, 
indented to the right of z = 0. Using the binomial expansion on (z — 2) 2n 

In 

Z - 



2) 2n = ( l) m (2n)\ ^2n-m z m 



m=0 



m\(2n — to)! 



and after some simplification, (14.31) leads to 

2e 2 ( — l) m (2n)! ^2n-m e (n+l)X 

' ^— ' m!(2n — to)! 



p(t\x) 



n=0 m=0 



1 

2tt7 



C" 



-dz. 



(4.6) 



where 



and 



2n + 1 

An = X + — 

2 2 



iV = 2ri — to + 2. 

We express the integral in (14. 6 j) in terms of parabolic cylinder functions, using 



2ttz 



— / ^e"^ +22/2 ^ 



1 



Br 



thus obtaining 

p(t\x) 



^7T 



n=0 



(n+i)x p -4/4 ( l) m (2n)! , , {2n _ m)/2 



m\(2n — to)! ^ 

m=0 v 7 



m—2n- 



l(z n ) + v29D m _2n-2(4) 



(4.7) 



where 



A n (2n + l)X + 6 



Replacing 2n - m by /, (g3D leads to (12T2H . 

If we let X — > and 0^0 with X/ y/Q (thus z n ) fixed, the term with 
I = in (I2.2ip dominates and we have 

2e 2 



p(t\x) 



' 7T 



e 



J2 e (n+l)X e -z^ Do{ 



2e 2 



n=0 

00 

E< 

n=0 
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-(2n+l) 2 X 2 /(40) 



(4.8) 



Here we note that X 2 /G = Z 2 JT and we used the fact that D (w) = e _u,2//4 . 
Since (14.81) is the same as (I2.18p . we have shown that case 5 is really a special 
case of case 6 in Theorem 2. 

Alternately, we can treat the problem on the (X, 0) scale by using the 
conformal map in section 2 to get (I3.15p . The poles of the integrand in (13.151) 
satisfy 

-2^/pX sinhr/' 



exp 



Since p = 1 — e, we must scale r\ to be 0(e). Let r\ = r\\e + r^e 2 + 0(e 3 ), we 
find that 



2t?i + 1 



-2X Vl 



(4.9) 



27/1-1, 

and r] 2 = Tji/2. Setting rji = u + iv we find that all the roots of (14.91) are on 
the imaginary axis, and with u = (14.91) becomes 

(Av 2 - 1) sin(Xv) = Avcos(Xv) 



or 



cos 



Xv 



2v sin 



Xv \ 



2v cos 



Xv 



+ sin 



Xv 
~2~ 



0. 



Thus we have either cot(Xt>/2) = 2v or tan(Xv/2) = —2v. Denoting the 
n th positive solution by v n = v n (X) we have v n (X) G ((n — l)ir/X, rnr/X), 
n > 1. The first solution satisfies cot(X^!/2) = 2t>! and has the asymptotic 
expansions 

71 AlT 167T 

x~x 2 + Jp 



MX) 



v x {X) 



1 



+ 0(X- 4 ), X -f oo, 



— + ^x 3 / 2 + 0(X 5 / 2 ), X -> 0. 
24 5760 V ; ' 



The higher order solutions behave as 

1 



V n (X) = U7T 



4 16 

X~ X 2 + X3 



+ 0(X~ 



X — > oo, 



and 



^n(X) 



7T 



71-1 — + 



X 



X (n-l)7r (n-l) 3 vr 3 



+ 0(X 2 ), ti>2,X^0. 
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Then in the s-plane the singularities satisfy s = 0(e 2 ) with 



r*(X;v n ) = (-1 - p + 2 v /pcosh?7)/e 2 

= -(^ + l/4)-(^ + l/4) e /2 + 0(e 2 ), (4.10) 

and we also have, at the singular points, 

exp [(1 - 2 v ^coshr ? + y/pe^) X/e] ~ e x/2 ~ iXvn . (4.11) 

Expression ( 12. 22ft follows by indenting the integration contour in ( 13.151) around 
all the singularities, using the symmetry (v — > —v) of the integrand and also 
using (I4.10p and (14.111) . Note that in heavy traffic all the singularities are 
close to 1] = and on the imaginary axis. If we retain only the leading term 
in r*, ( I2.22p is equivalent to (I2.20p . as can be seem by setting £ = r/f in (I2.20jl . 
with which the singularities satisfy (14. 9ft . 

If we consider even larger time scales, with t = a/e 3 = 0(e~ 3 ) (thus t = 
0(e~ 3 ) ), then the two singularities at rji = ±ivi dominate. Here V\ = V\(X) 
is the unique root in the interval (0,ir/X) of (12.231) . This leads to (12.24)) . 
Note that to compute the leading term on the a-time scale, we need the 0(e) 
correction term to r*. 

To compute the unconditional density p(t) on the scale t = 0(e~ 3 ) = cr/e 3 , 
we use ( 11.2)) with \i = 1 and ( I2.24H . For simplification, we write V for v±(X), 
and we have 

p(t) ~ e / F{X; V)e r * {x '' v)a/e e- x/e dX 
Jo 

poo 

~ e F{X-V)e^ y2+l l^/\\~ x - {v2+1, M^dX. (4.12) 
Jo 

This is a Laplace type integral. We let $(X) = -X - (V 2 + 1/4) er and 
must find the maximum of this function over X. Then $'(X) = — 1 — 
2aV(X)V'(X), and from (12.231) by implicit differentiation we obtain 

(1 + AV 2 )V 

V'{X) = 



[l + AV 2 )X + A 



The condition $'(X) = defines X = X(a) implicitly, where the inverse 
function is 

X(l + W 2 )+4 

° = 2V\l + W-) ' (413) 
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where V = V(X). Denoting the right-hand side of (I4.13P as k(X), we can 
verify that k'(X) > 0, so that k(X) is a monotonically increasing function, 
and that k(X) — > oo as X — > oo. From our asymptotic results V ~ \ j\[X 
as X — > 0, so that k(X) — > as X — ► 0. Hence there is a unique positive 
root of fc(X) = cr, which we denote by X, so that 

X(1 + AV 2 )+A ~ , 

a = -A^— 4^— , V = VYX). 4.14 

2\/ 2 (l + 4V 2 ) 

Then we use the standard Laplace method in (14. 12ft and get 



p {t) „ = F(X- y) e -(^+i/4W2 e *WA (4.15) 

-$"(X) 



The first root V\(X) = V satisfies 



cot ( — ) = 2V. (4.16) 



Letting ip = XV, we write ( 14. 16ft in parametric form, with X = 2%p tan(^/2) 
and V = | cot ("0/2). Then ( 14.141) becomes 

a = 4tan 3 (- ) (^ + sin^) (4.17) 



% 2, 

and we can rewrite (14.151) in terms of o and ip, using 



- $(X) = -X - (V + l/4)a = 2^ tan esc 2 ^| ) = F (V), (4.18) 



COS' 



3 (I) 



and 



F(X; V)e~ x/2 = V2/ , (4.20) 

sin (^J (-0 + sin^) 

$"Oh = 3 cot (|) (^ + sin^)+4cos 4 (f) 
1 J 2(^ + sin^) 2 ' 
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The expressions in fl4.17p - fl4.19l) are the same as those in Morrison [I] (corre- 
sponding respectively to (A.21), (A.16) and (A.17) in [4]). Using (Q7ll - (jOIj) 
in f)4.15p . we have 

^ ^ /2 ^(f) p _ FoW/£p _ FlW 

and this agrees with the result of Morrison (see (3.35) in [4J) if we note that 

p(t) = -— Pr[V > t) ~ ^(VOV'V) Pr[V > t). 

Finally, we discuss the connection between the unconditional tail expan- 
sion of Flatto (for fixed p < 1, see (1.3) in [8]) and the heavy traffic result of 
Morrison (see (3.35) in [1]). If p « 1 we have 

Pr [V PS > t] = i Pr [Wros > *] ~ Pr [W ROS > t] . (4.23) 

We expand the cr-scale result in jl] for a — > oo. It follows that ^ — > 7r as 
cr — > oo, with 

, / 7T — ?/A 327T . , „ 

a „ 4xa #\-JLy j — w . (4.24) 

Thus ^ ~ 7t - (327r) 1/3 <7 -1 / 3 . We replace i by pi in (3.35) in @j, so that 
a = e 3 pt. Then for a = e 3 pt — > oo the results in j3] become 

=3 'tt\2/3 



F (V)~ j/i* + 3(-J e(/it) 1 / 3 -4, 
^iW~^-i(|) 2/3 « 1/3 +2, 



and 

1/2 



2(sin^) (^y) - 2 14 /3 3 -V 27r 5/ 6e -3 (M) -5/ 6 _ 
Thus for cr — > oo the heavy traffic result in [1] becomes 



Pr|V > 4 ] ~ ^Llz^_2M , (4. 25 ) 
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where 

a* = 2 14 / 3 3- 1 /V/V 3 e 4 / e - 2 , 




If we let p = 1 — e and replace t according to At — (1 — e)\ii in Flatto's 
result ((1.3) in [8J), and expand that expansion for e — > 0, we also get (I4.25p . 
This shows that the interesting tail structure in (1.3) is also inherent in the 
results in [1], if we consider the latter for times even greater than the largest 
natural time scale. 



5 Discussion 

To summarize, we have derived several asymptotic formulas for the con- 
ditional sojourn time distribution in the M/M/l-PS model. Our analysis 
considers all of the space/time scales inherent to the problem, and leads 
to qualitative insights into the structure of the distribution. For example, if 
p < 1, x is large and t/x ~ 1/(1 — p), formula (12.51) simplifies to the Gaussian 



, I x (1"P) 3/2 

p{t\x) » - exp 

2^/pnx 



(1 — p) ( x 



t 



Apx \ 1 — p 



(5.1) 



which gives the spread about the well known mean value E [V(x)] = x/(l—p). 

In contrast, in the heavy traffic case, for x = 0(1) (short to moder- 
ate jobs) the sojourn time distribution is approximately exponentially dis- 
tributed on the large time scale T, as shown by (12.151) . We can show that 
this exponential distribution still holds for x = Zj^fi = 0(e -1 / 2 ) (where 
now the sojourn time must be scaled to be 0(e -3 / 2 )), but for very large 
jobs, where x = X/e = O^ 1 ), the distribution is concentrated on time 
scales t = 0/e 2 = 0(e -2 ), where part 6 of Theorem 2 applies. We also note 
that when the condition on x is removed, which corresponds to multiplying 
e -T/x e -x anc j integrating, we obtain the modified Bessel functions that 
appeared in [1]. 

Our results easily reproduce the unconditional asymptotics in [1] and [8] , 
and lead to a better understanding of what ranges of x are important for 
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removing the condition via (11.21) . We are presently investigating the more 
general M/G/l-PS model, and our preliminary results suggest that the basic 
scales for the conditional sojourn time asymptotics are very similar to those 
here, but that the unconditional asymptotics are quite sensitive to the tail 
behavior of the service distribution. 
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Table 1: The dominant singularity (rj — u + iv with v > 0) for various values 
of x and p < 1. 



X 

p 


0.01 


0.05 


0.1 


0.5 


1.0 


0.1 


6.5871 + in 


4.9639 + in 


4.2532 + in 


2.4866 + in 


1.5299 + ?'tt 


0.2 


5.8796 + in 


4.2448 + in 


3.5186 + m 


1.5625 + in 


2.3370? 


0.3 


5.3831 + i?r 


3.7333 + in 


2.9861 + in 


0.5195 + in 


1.6111? 


0.4 


4.9614 + in 


3.2909 + in 


2.5125 + in 


1.9978? 


1.2560? 


0.5 


4.5641 + in 


2.8626 + in 


2.0323 + in 


1.5312? 


1.0129? 


0.6 


4.1575 + in 


2.4049 + in 


1.4685 + in 


1.2118? 


0.8218? 


0.7 


3.7037 + in 


1.8498 + in 


0.4973 + in 


0.9534? 


0.6566? 


0.8 


3.1312 + i?r 


0.9559 + in 


1.8187? 


0.7188? 


0.5006? 


0.9 


2.1965 + in 


1.6812i 


1.1094? 


0.4748? 


0.3334? 


^•^^ X 

P 


2.0 


3.0 


4.0 


5.0 


10.0 


0.1 


2.2245?' 


1.5826i 


1.2849i 


1.0982?' 


0.6693?' 


0.2 


1.3682? 


1.0513? 


0.8736?' 


0.7553?' 


0.4694? 


0.3 


1.0394i 


0.8127i 


0.6810i 


0.5920?' 


0.3729?' 


0.4 


0.8359?' 


0.6599?' 


0.5562?' 


0.4855?' 


0.3099?' 


0.5 


0.6860?' 


0.5454?' 


0.4619?' 


0.4049? 


0.2620?' 


0.6 


0.5633?' 


0.4506?' 


0.3834?' 


0.3373?' 


0.2216?' 


0.7 


0.4544?' 


0.3655?' 


0.3123?' 


0.2759?' 


0.1842? 


0.8 


0.3492?' 


0.2823?' 


0.2423?' 


0.2149?' 


0.1461?' 


0.9 


0.2342?' 


0.1903?' 


0.1641? 


0.1461?' 


0.1012?' 
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Figure 1: The contour C. 
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Real_part 
Imaginarypart 

Figure 2: The singularities of the integrand in (13.151) for x — 1, p = 0.3. 
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Real_part 

Imaginary_part 

Figure 3: The singularities of the integrand in ( 13.151) for x = 0.5, p = 0.3. 
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Sample text for slashbox.sty 

The usage is pretty straightforward, such as 



Date 

Room 


5/31 


6/1 


6/2 


6/3 


6/4 


Meeting Room 












Auditorium 












Seminar Room 













You may include a newline (\\) in 'Room' and/or 'Date'. Note that you 
will get spaces aside the slash line if there is a wider column in the same 
column of a different line. In such a case, you need to specify the width of 
the slashed column by saying 



Date 

Room 


5/31 


6/1 


6/2 


6/3 


Room With a Long Name 










Auditorium 










Seminar Room 











The specified width is neglected if it is narrower than the natural width 
of the column. 

\ (back) slashbox assumes by default that there is a blank space of width 
\tabcolsep on both sides of the column. Thus the slash line might exceeds 
the boundary when you use @0 etc. 

You can avoid it by specifying 



Date 
Room ^^^-^^ 


5/31 


6/1 


6/2 


6/3 


6/4 


• Meeting Room 












• Auditorium 












• Seminar Room 













Here [1] tells the command that there is no extra space on the left of 
this column. You can use [r] and [lr] likewise. You have to also specify 
the width of the column in this case, but it can be Opt. 



